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Abstract 

We study the Weyl representation of metaplectic operators suggested 
by earlier work of Mehlig and Wilkinson. We give precise calculations 
for the associated Maslov indices; these intervene in a crucial way in the 
Gutzwiller formula of semiclassical mechanics. 



1 Introduction 

In an interesting paper |Hj the physicists Mehlig and Wilkinson introduce, in 
connection with their study of the Gutzwiller semiclassical trace formula, a 
class of unitary operators S : L^(R") — > L^(K"). These operators are defined 
as follows: let S G Sp{n) have no eigenvalue equal to one; to S one associates 
the Weyl operator 



2ttJ y^\det{S-I) 

where T{zo) is the Weyl-Heisenberg operator and 

Ms = y{S + I){S - I)-^ (2) 

/ being the identity and J the standard symplectic matrix (see below). The 
index v is an integer related to the sign of dot (5* — /), and which is not studied 
in the general case in Mehlig and Wilkinson moreover show that 

R{SS') = ±R{S)R{S') (3) 

for all S, S' for which both sides are defined. They claim that these operators 
belong to the metaplectic group. This property is however not quite obvious; 
what is acceptably "obvious" is that R{S) is a multiple by a scalar factor of 
modulus one of any of the two metaplectic operators ±5* associated to Mp{n); 
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this is achieved using the metaplectic covariance of the Heisenberg-Weyl opera- 
tors (see below) . The purpose of this paper is to precise Mehlig and Wilkinson's 
statement by comparing explicitly the integer j/ in J^l with the Maslov indices 
on the metaplectic group we have studied in a previous work |5] . This is indeed 
important -and not just an academic exercise- since the ultimate goal in |S] 
is to apply formula Q to give a new proof of Gutzwiller's trace formula for 
chaotic systems. It is well-known that the calculation of the associated "Maslov 
indices" is notoriously difficult: it suffices to have a look on the impressive bibli- 
ography devoted to that embarrassingly subtle topic. We will, in addition, give 
a semiclassical interpretation of R{S), expressed in terms of the phase space 
wavefunctions we introduced in [2 0] . 

Remark 1 An alternative approach to the results of this paper would be to use 
Howe's beautiful "oscillator group" method '5/ (see 'J] for a review); this would 
however in our case lead to unnecessary technical complications. 

Notations 

We denote by a the canonical symplectic form on M^" = M" x 

tT(z,z') = {p,x') — {p',x) if z = {x,p), z' — [x'p') 

that is 

cr(z,z') = (Jz,z') , J = J Q . 

The real symplectic group Sp{n) consists of all linear automorphisms S : R^" — > 
R^" such that a{Sz, Sz') — cr{z, z') for all z, z' . It is a connected Lie group. We 
denote by €x and ip the Lagrangian planes R" x and x R^, respectively. 
iS(R") is the Schwartz space of rapidly decreasing functions on R", and its dual 
<S'(R") the space of tempered distributions. 



2 Prerequisites 

2.1 Standard theory of Mp{n): Review 

The material of this first subsection is quite classical; see for instance 013] and 
the references therein. 

Every 5* £ Mp{n) is the product of two "quadratic Fourier transforms" , 
which are operators Sw,m defined on S{X) by 

Sw,mf{x) - (^^) J e^'^(-'-')/(x')d"a:' (4) 

where Vl^ is a quadratic form in the variables x,x' of the type 

W{x,x') = ^{Px,x) - {Lx,x') + ^{Qx',x') (5) 
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with P = , Q = , detL ^ 0. The integer m appearing in corresponds 
to a choice of argdetL: 

mTT = arg det L mod 27r 

and to every W there thus corresponds two different choices of m modulo 4: if 
m is one choice, then m + 2 is the other (this of course reflects the fact that 
Mp{n) is a two- fold covering of Sp{n)). The projection tt : Mp{n) — > Sp{n) is 
entirely specified by the datum of each TT{Sw,m), and we have TriSw,m) = Sw 
where 



{x,p) = Swix',p') 
In particular, 



p = d.j:W{x,x') and p' ^ —dxiW{x,x'). 



PL-^Q - 



(6) 



is the free symplectic automorphism generated by the quadratic form W; observe 
that Sw^p r\ £p = for every W. The inverse — S^.^ of Siy^m is the 

operator Sw',m' where W*{x,x') — —W{x',x) and m* — n ~ m, mod 4. Note 
that if conversely S" is a free symplectic matrix 



S ■ 



e Sp{n) , detB^O 



(7) 



then S = Sw with P = p-^A, L = B-\ Q = DP- 



2.2 Heisenberg— Weyl operators 

For Zq = (xq,po) we denote by T(zq) the translation z i — > z + zq; it acts on 
functions by push-forward: T(zo)/(z) — f{z — zq). We denote by T(zo) the 
corresponding Heisenberg- Weyl operator: for / e we have 

f (zo) = e^''P"'^^-^^P"''="^^f{x - xo). 

The operators T(zq) satisfy the metaplectic covariance formula: 

Sf{z) = f{Sz)S (S = 7r{S)) (8) 

for every S € Mp{n) and z. In fact, the metaplectic operators are the only 
unitary operators, up to a an factor in satisfying (|SJ|: 

For every S € Sp{n) there exists a unitary transformation U in 
satisfying ^ and U is uniquely determined apart from a 
constant factor of modulus one. 

The Heisenberg- Weyl operators moreover satisfy the relations 

f (zo)f (zi) = e-^'^(^°^^^)f (zi)f (zo) (9) 

f (zo + zi) = e-^(^«-^^)f (zo)f (zi) (10) 
as is easily seen from the definition of these operators. 
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2.3 Weyl operators 

Let a™ be the Weyl operator with symbol a: 

= (27)"/ e^^^^--y^a(\{x + ylp)f{y)crycrp; 
where / £ 5(R"')equivalently 



where is the symplectic Fourier transform F„a defined by 
The kernel of a™ is related to a by the formula 



The Mehlig- Wilkinson operator is the Weyl operator with twisted Weyl 
symbol 

a{z]=(—\ (Ms^o.zo) QlN 

2.4 Generalized Presnel Formula 

We will use the following formula, generalizing the usual Fresnel integral to 
complex Gaussians. Let M be a real symmetric nxn matrix. If M is invertible 
then the Fourier transform of the exponential exp(i (M x, a;) /2) is given by the 
formula 

where sgnM, the "signature" of M, is the number of > eigenvalues of M 
minus the number of < eigenvalues. 

For a proof see for instance 1;, App. A. 



3 Discussion of the Mehlig— Wilkinson Formula 

The Mehlig-Wilkinson operators R{S) are Weyl operators with twisted Weyl 
symbol 

We begin by giving two straightforward alternative formulations of these oper- 
ators. 
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3.1 Equivalent formulations 

We begin by remarking that the matrix Mg = ^J{S + T){S — I)~^ is symmetric; 
this immediately follows from the conditions 

S e Sp{n) <=^ S^JS = J ^ SJS^ = J. 

Notice that lO can be "solved" in S, yielding S = (2M — J) ^(2M + J). 

Proposition 2 The operator 

R{S) ^ (—] , / e5<A'^s^"'^")f (zo)d'"zo (13) 

can be written in the following alternative two forms: 

R{S) = (^^^ zV|det(5-/) J f (5zo)T(-zo)d'"^o (15) 

for det{S -1)^0. 
Proof. We have 

i J(5 + I){S - /)-! - i J + J(5 - 

hence, in view of the antisymmetry of J, 

(Mszo, 2o) = {J{S - ly'zo, zo) = a((5 - ly'zo, zq) 

Performing the change of variables zq i — > (5 — /)~^zo we can rewrite the integral 
in the right hand side of (|13|l as 

g|(Ms^o,^o)f(^)^2„^^ = / e5'^(^''^(^--^)^'')f ((S* - /)zo)rf^"zo 

hence dJ. Taking into account the relation (|10|) we have 

f ((5 - I)zo) = e-5-(s^'"^'')f (5zo)r(-zo) 
and formula l|15(l follows. ■ 

Corollary 3 We have R{S) = csSwm where c is a complex constant with 
|c|-l. 

Proof. We begin by noting that R{S) satisfies the metaplectic covariance rela- 
tion 

R{S)f{zo) = f{Sza)R{S) 

as immediately follows from the alternative form (|15|l of R{S). On the other 
hand, a straightforward calculation using formula (|14() shows that R{S) is uni- 
tary, hence the claim. ■ 
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3.2 The case S = 



W,m 



'A- I B ' 




B 




C D-I 




I D-I 





We are going to show that the Mehlig- Wilkinson operators coincide with the 
metaplectic operators Sw,m when S = Sw and we will thereafter determine the 
correct choice for v] we will see thast it is related by a simple formula to the 
usual Maslov index as defined in 

Let us first prove the following technical result: 

Lemma 4 Let Sw be a free symplectic matrix We have 

detiSw -I) = det B Aet{B-^ A + DB'^ - - (B^)"^ (16) 

that is, when S is written in the form 0); 

detiSw -I)= det(i-i) det(P + Q-L~ lF). (17) 

Proof. We begin by noting that since B is invertible we can write S — I as 

\C-{D-I)B-\A-I) 6 
B-^{A-I) I 

hence 

Aei{Sw -I) ^ det B detiC - {D - I) B-\A - I)). 

Since S is symplectic we have C — DB^^A = —{B'^)^^ (use for instance the 
fact that S'^JS = SJS'^ = J) and hence 

C-{D- I)B-\A - /)) = B-^A + DB-^ - B'^ - (5^)"^ 

the Lemma follows. ■ 

Proposition 5 Let S be a free symplectic matrix ^ and R{S) the correspond- 
ing Mehlig-Wilkinson operator. We have R(S) = Sw,m provided that v is chosen 
so that 

v = m~lnert{P + Q-L-L'^) mod4 (18) 

(Inert (P + Q — L — L'^) the number of < eigenvalues of the symmetric matrix 
P + Q-L-L^). 

Proof. Recall that we have shown that R{S) — csSw.m where cs is a complex 
constant with \cs\ = 1. Let us determine that constant. Let 5 £ 5'(R") be the 
Dirac distribution centered at a; = 0; setting 

1 " 



^\dei{Sw-I) 



we have, by definition of R{S) 
R{S)S{x) = C 
= C 



ei{Ms{x,po),{x,Po))^i{p,x)gl^^ xo)d^'''ZQ 
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hence, setting x — 0, 

R{S)S{0) = C I e*<*^^(°'^"')'(°'f«)><5(-xo)d2"zo 



that is, since / S{—xo)d^xo = 1, 

Let us calculate the scalar product 

(M5(0,po), (0,Po)) - ^((^ - /)-'0,po), (0,Po)). 

The relation = (5 — 7)^^(0, po) is equivalent to S{x,p) = {x,p + po) that 

is to 

p + Po = dxW{x,x) and p ~ —dx'W{x,x). 
Using the explicit form Q of together with Lemma 0] these relations yield 

x={P + Q~L~ L^)-^po ; p={L-Q){P + Q-L- L^y^po 

and hence 

(Ms(0,po), (0,po)) = -{{P + Q-L- L^)-'po,po) . (20) 
Applying Fresnel's formula (|12|l we get 



since 

- I det L| V2| det(P + Q - L - L^jr^^ 

V|det(5-/) 
in view of H17|l in Lemma 01 we thus have 

R{S)d{0) = (^-Ly\-e-T-gn(^+Q-i-i")|detL|i/2. 

Now, by definition of Sw,m, we have 

Sw,m6i0)^ (^^) z™""/2|detL|i/2 



/2. 



i'^-^ sgn(P+Q-L-L-' ) _ •m-n/2 



hence 

It follows that we have 

1 n 
1/ — - sgn(P + Q — L — L ) = m — — mod 4 

which is the same thing as (|18|l since P + Q — L — has rank 
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3.3 The general case 

Recall that we established in Lemma 0] the equality 

det{Sw -I)^ det det(P + Q-L-LF). (21) 

valid for all free matrices Sw G Sp{n). Also recall that every S G Mp{n) can 
be written (in infinitely many ways) as a product S — Sw,mSw' ,m' ■ We are 
going to show that Sw,m. and Sw',m' in addition always can be chosen such 
that det{Sw,m — I) ¥^ ^ and de\.[Sw',m' — I) 7^ 0- For that purpose we need the 
following straightforward factorization result (see [3j): 

Lemma 6 Let W be given by |3|); then 

Sw.rn - V-pML,mJV^Q (22) 

where 

V-Pfix) = ei(P-'-^fix) ; ML,,nf{x) = *'"x^^^/(Lx) 
and J is the modified Fourier transform given by 




Let us now state and prove the main result of this section: 

Proposition 7 Every S G Mp{n) is the product of two Mehlig-Wilkinson op- 
erators; these operators thus generate Mp{n). 

Proof. Let us write S = Sw,mSw' ,m' and apply H22|l to each of the factors; 
this yields (with obvious notations) 

S = V-pML,,nJV_f^P-+Q)ML',^'JV^Q'. (23) 

We claim that Sw.rn and Sw',m' can be chosen in such a way that (let[Sw,m — 
I)^Q and Aci{Sw' ^m- - I) 7^ that is, 

det(P + Q-L-L^)^0 and det(P' + Q' - L' - L'^) 7^ 0. 

This will prove the assertion in view of H21|l . We first remark that the right 
hand-side of (|23|l obviously does not change if we replace P' by P' + \I and 
Q hy Q — \I where A G M. Choose now A such that it is not an eigenvalue of 
P + Q — L — and —A is not an eigenvalue of P' + Q' — L' — L'^; then 

det(P + Q - \I - L - L^) and det(P' + \I + Q' - L - L^) ^ 0. 
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